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ABSTRACT 



Submitted to the Department of Naval Architectiu’e 
and Marine Engineering on May 17 > 1968 in partial ful- 
fillment of the requirements for the Master of Science 
Degree in Electrical En,gineering and the Professional 
Degree, Naval Engineer. 

This paper presents a procedure for the design of 
superconducting electrical machines. The magneto- 
static field problem of a cylindrical superconducting 
machine is solved. Equations for the magnetic fields 
are developed, and from these, expressions for the 
various machine parameters are obtained. These expres- 
sions arc adapted to a computer solution of the design 
problem. The computer program obtains the design 
parameters for a minimum volume machine design. 

A design study for a typical marine electrical 
propulsion system is conducted. The results of this 
study indicate quantitatively the reduction in ueight 
and space to be obtained in using superconducting 
electrical machines vice conventional electrical 
machinery. The study also presents some general 
characteristics of superconducting machines . 
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I. IKTRODUCTIOW 



A. BackRroi.md 

The relative flexibility and simplicity accruing to 
electrical power systems, in comparison to mechanical power 
systems, has long been recognized. Specifically, electrical 
power, and transmission systems provide flexibility of 
installation, are easily suited to automation and have a high 
degree of dependability. 

Despite these advantages, and others that could be obtained, 
the full benefit of electrical povrer systems has not been realized 
in marine propulsion applications. To be sure, electric drive 
systems have been tised in ship propulsion, however, such use has 
been mostly limited to small or intermediate size propulsion 
plants. The primary reason that electric drive systems have not 
been adopted for large power plants has been that the weight and 
space for such systems is greater than that required for geared 
propulsion systems. In addition, increased cost (as much as 2Qffo) 
and lox'ier efficiency (85^ for DC) have been contributory factora^^^ 

It is advantageous to consider some of the benefits that 
could be obtained vrith electric propulsion in order that aspects 
of its utilization might be judged in proper perspective. 

As mentioned above, an electric propulsion system affords 
greater flexibility of installation. In constrast with a 
conventional geared propulsion system, extensive shafting is not 

"riuvibers in parant}icses indicate applicable refex'cnce in the 
b ib 1 5 og r aph y . 
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required and the location of the pov:or sooipce is not as 

restricted. In addition, the large, noisy reduction gear could 

be eliminated. The noise generated by this component is of 

increasing military significance. 

If a cycloconvertor" is used, additional benefits derive. 

It would be possible to use a constant frequency generator, 

thus permitting a more efficient design of the prime-mover. 

This will also permit greater use of nevrly developed povrer 

sources such as the gas turbine. Added to this would be more 

adaptation of remote and automatic controls, plus much improved 

(2 * 5 ) 

maneuvering characteristics.' 

Thus it seems that if the disadvantages (weight, space, and 
possibly efficiency) of electric propulsion are reduced, it 
vjill be given greater consideration as a:i altei’native to geared 
propulsion in larger pov/or systems. 



® * The Role of Superconductors 

It has been realized for some time that when the temperature 
of an electrical conductor is reduced to a fevr degrees absolute, 
the rosistanoo of the conductor decreases to an immeasurable 
amoimt. The implications of this fact have considero.ble value 
for uso in electrical machines. If the resistance in the lunding 
of an electrical generator or motor is removed, much greater 
cui’rent densities can be achieved with an accompanying increase 



in the magnetic flx’sc densities within the machine. 



( 3 ) 



If this 



"x'n electronic povjoi'’ tranem:-,ssion device i 
voltage and frequency. See reforoneo Ii.. 



-ha,t affords changes in 



- 3 - 



increase in flixx density is obtained, then greater power can 

be generated within a given volume. Thus it follovrs that the 

volume of an electric machine can be reduced if superconducting 

( 1 ) 

windings are used. Studies to date, have demonstrated the 
significant reduction in volume and vzeight in superconducting 
machines , 

However, the application of superconductivity is not vrithout 
its complications. The first question vrhich must bo asked 
concerns the method of operating a machine vrhcrein all, or at 
least some, of its parts are at a temperature in the region of 
liquid heliTAm. 

It is at this point convenient to limit the scope of the 
examination. The ensuing discussion will treat only machines 
with a rotating field trinding, operating in the superconducting 
region, enclosed by a stationary armature at ambient temperature. 
This particular configuration is only one of several that could 
be considered. 

There are two particular reasons why this configuration is 
selected. Tlie more obvious is that it vjoxild bo irapractical to 
transfer the high machine power output through slip rings, as vrould 
bo necessary if the armature were rotating. The other reason 
results from the fact that large heat losses are cncoimtered 
when superconductors ai’e subjected co alternating currents, 

Thus, only the field vrinding is made superconducting. 



Nonetholcss, most of the techniques in this paper are 
generally applicablo to any cylindj'ical elocbrical device. 



I 
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With the discussion constrained to considering a rotating 
superconducting field winding, it is possible to dilineate some 
design procedure. Thus far, it has been found that the \iso of 
a thermal bottle, or dexirar", filled xxith liquid helium and 
encompassing the rotor, is a satisfactory solution to the 
problem of maintaining the vrindings in a superconducting state. 
Several techniques are under development for introducing the 
helium into the rotating dewar and passing electrical leads 
into the low- temperature region. Generally it is favorable 
to use the liquid helium dexxar in conjunction x-;ith a liquid 
nitrogen dexirar irmnediately exterior. The helixmi maintains the 
rotor at about 1}..2°K x^hile the nitrogen offers a buffer 

for improved cryogenic performance. At first, one might think 
that the refrigeration load for maintaining such low temperatures 
vjould be excessive. Ilowover, it should be realizod that once 
the lovT temperatxu'es are obtained, there is very little heat 
generated by the windings. Indeed, heat introduced into the 
low- temperature region tlirough the shaft and support structure 
is more significant. A reasonable estimate of the pov;or 
requirements for refr3.gcration during operation of a super- 
conducting machine would bo about 0,5 percent of rated power. 

In addition to the requiromont that a superconductor be 
maintained at a lex; temperature, there is a limitation on the 
amount of current carried v;hile exposed to a magnetic field. 

As indicated in Figure 2, if the magnetic flux density increases, 



'Ikuned after the inventor. 
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the permissible current density is reduced. If current or 
field exceed the limits imposed (the so-called critical 
values), the conductor reverts to its normal state. Thus 
it is necessary, in any particular machine design, to 
determine that combination of maximum flux density and 
current density which best suit the design. 




Figure 2 Current -Flux Relationship for a 

Typical Superconductor 

The very fact that a s\iperconduct ing machine develops 
extremely intense magnetic fields invites yet another 
complication. Since fields of the order of several kilo- 
gauss can be expected vrithin the machine, disturbance of 
equipment in the vicinity of the machine will certainly 
result, as vrell as unbalanced loads on the superconducting 
field wiiiding. Therefore, it is necessary that a shield be 
provided to confine th.e magnetic fields vrithin the machine. 
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This shield would typically be of laminated steel con- 
struction and x>;ill be a significant portion of the sizo 
and w’eight of the machine. 

C . Design Program 

Though superconducting machines of the size suggested 
by the foregoing discussion have not yet been built, they 
are certainly feasible. Anticipating the construction of 
such machines, it is considered desirable to investigate 
some of their expected properties. Since there is vested 
interest in the reduced size and X'xeight of these machines, 
it is natural that a design procedxare be developed for 
obtaining minimum weiglit and volurae. 

The objective herein is to find a procedure for finding 
design characteristics of least size machines. The procedure 
developed is founded on a field-theory sol.ution of the 
electric machine and is adapted to a computer solution. 
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II . PROCEDITP.S 



A, Developinoiit of a Model 

The first step in the analysis was to obtain a 
mathematical model for the magnetic fields involved. This 
model i-7as necessary to determine the parameters vjhich would 
be used in machine design. As mentioned previously, the model 
and attendant equations are for a rotating field winding inside 
a three-phase stator winding. 

Analysis proceeded by considering both the field and stator 
windings to consist of cylindrical current sheets. Equations 
for the magnetic fields generated by each vrlnding vrero developed 
separately and later combined to obtain necessary inductance and 
energy expressions. 

To illustrate the analysis, consider the field x/indings’. 

As illustrated in Figure a cylindrical current sheet, K(r). 
of thickness dp , located at some radius r, R,ZrZR^, is examined 
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Fourier analysis of this current sheet yields: 



K(r) 



E 



n odd 



M^dr 

"nir 



sin np 0 , 



(II-l) 



Knovring that these current sheets exist betxireen r ” and 
r = and that the fields associated x:ith each sheet are 

( 

Laplacian, enables determination of the field expressions, 
After finding an expression for the magnetic field for r ~ R 2 » 
one can find the contribution due to the shield. Since at 
r “ R , x^ith an infinitely permeable shield. 



II 



s0 




(II-2) 



and at r “ 0, 



H 



finite 






an expression is found for the fields induced by the shield. 
The expression thus obtaixiod is then added to that for the 
fields generated by the field x;inding. 

This same procoduna obtains for deterriini''.g t}ie magnetic 
field systoia resulting from the armat'aro vandingc , The only 
difference involved is in the current sheet model used, since 
the armature is r three-phase xrinding, X'.m.th 60-degroc phase 



belts . 
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Fourior analysis of the armature ciJirrent sheet yields 



K (r) 



a 



E itJ dr 



ntf 



n odd 



cos -j- sin np0 . 



(II-4) 



Expressions for magnetic field systems of the field and 
armature windings are presented in Tables 1 and 2. 

Using the equations for the magnetic fields ^ an expression 
for energy stored in the armature field system is obtained frora 



/'2rr 



Vf = 
a 








V-% ) IP tip d0 . 



J 



o 



o 



(I I- 5) 
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Equating the result in 





L 

a a* 



( 11 - 6 ) 



where 



I- = 



tt(R^ - 
o 

Til 



r2) 



(11-7) 



yields an expression for the armature self-inductance. A 
similar calculation is performed to obtain the field winding 
self- inductance . 

Mutual inductance is foimd from 



^f^af 




(II-8) 



v/here X is the total field-winding flux linlcage of the armatune 
windings. By recognizing that 6K /rr(R - R") expi'esses the 
number of atmature turns per unit cross-section, then the 
second spatial differential of the flinc linlcage is expressed 
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Kacrnctic Field FTuations for Field ’./indin; 
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TABL 5 2 

Ma,o:.netle FIs Id Souatloiis for 
Phase A of Armature V<'lndlnp;s 
np 2 
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integrating this expression over the cross-section of the 
armature yields the total flux linkage and thus the mutual 
inductance. Expressions for field and armature self- 
inductance and mutual inductance are presented in Table 3. 

Design Equations 

If one considers the araplitudes of the line-to-neutral 
armature voltage and armature current, one obtains 



= CO^KIj, 


(II-IO) 




(II-ll) 


' Sn » 




a 




= |v I , 

2 a a 


(11-12) 



where 

I r= "^^2 ” ^l^*^f (11-13) 

^ 2N^ 

This power is internally generated power, which must 
be reduced, to account for armature impedance and the load 
poxirer factor, in order to obtain rated output power. 
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From the vector diagram in Figure where 

( 11 - 14 ) 

the ratio V./V is obtained. Then rated povrer is 

"D 3. 

p = Pg(v^/v^) cos e. , (II-15) 



In developing the equations desci’ibing the machine, seme 
approximation vras necessary to account for the effects of the 
end-turns . Figure 6 shov:s the geometry involved in making 




If ono considers an outer 
a given phase vnnding, it 



armature turn at the mean radius of 
is soon thau the anguiai displacement 
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of this coil is 4.Tr/3p radians. Since this coil is at a radius 

r |(R^ + R^) (11-16) 

the approximation can be made that 

Al= I • ^ ■ |(R„. R.) = ?-{R„+ R^). (11-17) 

Using this expression, it has been found^^^ that a reasonable 
length to assiime for self-inductance calculations is the total 
length 

1^ = 1 + 2 Al, (II-18) 

whereas the length to be used for mutual inductance is 

1 ^ 1 -s- Al. (11-19) 

m 

^ ♦ Design Limitations 

As described in the Introduction, the lo^^; resistance 
character of a superconducting coil is not only a function of 
temperature, but also depends on its current density and the 
magnetic field to vjhich it is exposed. Thus for a given design 
it is necessary to determine the maxiircr magnetic flxu: density 
associated with any particu3ar current density, to ensure that 
the superconducting field winding v;ill not revert to its normal 
state. The particular mechanics employed to observe this 
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liBiitation are detailed in the program construction techniques 
of Appendix C . 

As also mentioned in the Introduction, the machine design 
considered employs a shield (constructed of laminated steel) 
which contains the dense magnetic fields iiiithin the machine. 
For this shield to be effective, the fields it constrains must 
be less than those v;hich v;ill cause saturation of the shield. 
If the shield material saturates at some density then 

one limitation is 



B , (r 
max. ' 



R ) - B . 
o rated 



( 11 - 20 ) 



Another limitation exists in the requirepient that the maximum 

flxix density carried by the shield must never exceed B . 

rated 

As illustrated in 




F igure 7 Shield Flux Pattern 
Figure 7, the total flux accumulated by the shield 
of n'/2p must not be siuTicient to cause satui’ation 



over an angle 
in the cross- 
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section at A-A. Since 



ID 



A-A 



/rr/2p 

o 






( 11 - 21 ) 



a final constraint is obtained lii 



*A-A - ®rated <"os ‘ %'>\' 



(II-22) 



There are two final considerations which nay govern a 

machine design. First, the heat dissipation in the armature 

v;inding must be considered. Generally the maJKimum permissible 

current will be determined from the limitations imposed by the 

method adopted for cooling the aimaturo. For purposes of this 

design it X'fas assumed that an average armature current density 
6 2 

of 1.6 X 10 amp/m'^ could be accepted with chilled v;ater cooling 
of the winding. Finally, one must consider the speed of the 
machine rotating parts. The strength of the rotor structure 
will provide a limit to the rotor tip speed that can be sustained, 
On consideration of previous conventional machine design it was 
determined that '(00 ft/sec xira.s a reasonable upper va3.ue for tip 
speed. Hov;ever, since methods of constr-acting superconduct5.ng 
vdndings are fa.r from being fixed, this value for maximum tip 
speed is speculative. The tip speed of the rotor t/ill also bo 
governed by the performance of the rotating dcvrar. In this 
respect, the fluidic and thermodynamic properties of liquid 
heliixm, itndsr high centrifugal forces, will govern, f-inco the 
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aspects of this consideration are yet to be detei’rained, it is 
assumed here that the previously stated raaximvim tip speed is 
acceptable. 



. 



I 
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III. RESULTS 



Tho significance of this v.»ork lies in the results of a 
design study which was conducted, employing the design 
procedure described in Part II. The design study involved 
what was considered a typical marine electrical propulsion 
system; a turbino driven generator rated at 26 MV/ (35^000 hp) 
at 3600 rpm, and two direct drive propulsion motors rated at 
17^500 bhp (13 MW) at 300 rpm. Other paraiaoters assuraed v;ore; 



■"f = 


D 0 

1.0 X 10° amp/m'^ 


J 

a 


P 

1.6 X 10 0.mp/m 


B " 

max 


2 

5.0 web/ra 


^rated 


2,0 v;eb/m^ 



Pj. - 40OO kg/m^ 

p = i{.000 kg/m^ 

^ (copper vrith .I|5 packing factor) 

p ~ 8000 kg/m^ 

(steol) 

S . 05m 



Tip speeds of up to 21I4. m/sec (700 ft/sec) wore considered 
in both the generator and motor designs. Also, for tho 
generator, designs v;ith po;isr factors of 1.0, .850 and ,70? 
vforo examined. Design characteristics of machines with 
1 to 6 polo pairs were comp'nted for all variati vons . 



# 

} 



i 



I 
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The computer results and graphs of the characteristics 
of some of the more interesting designs are included in 
the following pages. 

One quantity, not calculated by the computer pi'ogram, 
that should be examined is the machine efficiency. 

Since 



R “ Pr^^a^t 
a ; 



(III-l) 



where = resistivity = 1.9 x 10” ohm-m and cross-section 



area is 



, „ Tf(l-x^)R^ * .45 



(III-2) 



6N 



a 



for a .45 packing factor, then for the optimura generator 
design (pf = .70?^ P ” 6, = 2l4 m/sec) 

R /N^ - 6.7 X 10~^ ohjn/turn^. 
a a 



Knov/ing K I , it is possible to find 

S. Q. 

watts . 

These losses, ho^^rever, are for only one phase and vrere 

computed vrith a peak valvie for I instead of nns . Therefore, 

& 

resistive losses for all three phases are 
ifn ■' 965 * 4 1450 watts. 

If an equal araount is allowed for eddy current losses, a 
figure of 2900 watts is obtained for the total losses. 

Since this is for a 26 lit' generator, this 
a negligible reduction in efficiency. 



loss -will cause 



05'J V JF = U.lOE 09 A/H2 OMEGA = 377. RAD/SEC 
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Figure 8 

Generator Specific Weight 
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Figurg 9 

Motor Specific Weight 
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IV. DISCUSSION OF RESULTS 



The most significant result cf this study lies in the 
low values of specific weight (Ib/hp) for superconducting 
machines. As indicated in Figure 8 specific weights on the 



order of 0.25 Ih/hp" are theoretically possible in machines 
of the size typically found in marine propulsion plants. 

The results of the study of direct-drive propulsion 
motors i-fere less conclusive. It is apparent from Figure 9 
that specific v;ei{^ts of l.l-i-O Ib/hp are possible and that 
much lovjer specific weights are likely to be achieved in 
machines v;ith larger numbers of poles and higher tip speeds. 
Nonetheless a specific weight of I .40 is still an improvement 
over a conventional machine. This is evident from a 
comparison with the data for the propulsion motor in the 
USS Hunlcy, AS-31j given in Table 4 .. 



hp 



15x000 



pole pairs 



26 



rpm 



167 



Ib/hp 



10.9 
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be a limit in machine design. Hov;evor*, this is not likely 
to be the case in lovror speed machines such as the direct 
drive motor. 

Machine characteristics such as low and high efficiency 
are expected as a result of the smaller amount of ai^mature 
winding required. 

Aside from the significance of the data given for the 
various machine designs, there ^^rero two other notable 
results of this work. As described in the procedure for 
determining the field I'jinding design (detailed in Appendix B), 
in order to reduce the number of variables involved it via.3 
assumed that the fields induced by the shield did not affect 
B However, the study conducted in Appendix E indicates 

that in some instances (especially for small values of p) 
this assumption is not warranted. Prom Pigujae 16 it is 
found that for the optimum generator design (x ^ .822, 
y = .8Ij.8, p = 6) the error introduced is about one percent, 
vrhereas for the optimum motor (x , 801 , y = ." 887 , p = 6) 
the error is about .80 percent. These accuracies are ex- 
pected to be much greater than the accuracy of the math- 
ematical model being used. 

Aiiother result that becarde evident during the design 
study was that the linear interpolation teclinique used to 
find the min.1mum volume design (described in Appendix D) 
x/as not sufficiently accurate in all instances. Since the 
minima of the volume vs. R curve is ouJte broad for 







. 



I 






p = 1 and 2f the linear interpolation yielded erroneous 
results in a few designs. 

The accuracy of solution provided by the design program 
(within the limits imposed by the assvimptions as discussed 
in Appendices B and E) is arbitrary. In the design study, 
a variation of 1% in output power was allowed. Likewise, 
in instances of solving for maximum fl\ix density, the rate 
of change was required to be less than 10”^ at the calculated 
maximum. By changing these tolerences, the designer can 
obtain the accuracy desired in the solution. 
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V. CONCLUSIONS AND REG Oi-H-ffiNDAT IONS 

Prom the results of this study, the following can be 
concluded: 

1. Superconducting electrical machines should provide 
considerable savings of v/eight and space in electrical 
propulsion systems. 

2. The electrical efficiency of superconducting machines 
is substantially higgler than the efficiency of conventional 
machines . 

3. The solvition of the design problem for a super- 
conducting machine can bo made at least as accurate as the 
design of a conventional machine. 

ij-. High speed machines ^•;ill be limited in speed of 
rotation by the maximum tip speed that can be sustained by 
the rotating paj?ts. 

5. An optimum design for a slov; speed machi.ne will 
require a largo ntmibcr of poles. 

In the conduct of further investigation of the optimal 
design problem, or in application of the design program 
presented here, the follotiing is recommended: 

1. The body of data (array B(p, y) ) representing the 
graph of B . /const, vs. y, should be expanded to account 
for larger valiies of p. 

2. The design program should be modified to allov; 
variation in VTIP. As indicated by Figures 8 and 9j aji 
opt! um machine design may occur for some value 



of VTIP 
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that is less than the maxirnuin VTIP, By searching over 

a range of this variable, a nmnerical technique could 

easily bo used to find the optimum design. 

3, I'/hen a design is obtained, it is necessary to 

ensure that the error in B_ , as described in 

p max 

Figures 13 to 16, is not excessive. It may be desirable 
to incorporate these figures in the design program, as 
was done with array B(p, y). It vrould then be possible 
to check each design v.dthin the program, making appropriate 
adjustment in the design parameters until a satisfactory 
result is obtained. 

1|. To ensure greater acciaracy in the subroutine MIN, 

( 7 ) 

described in Appendix D, a false position solution teclmiquio ' 
should be sxibstituted for the linear interpolation miethod. 

A basic routine, such as the one presented on the follovring 
page, might be considered. 
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IF(STAR-1)11»12?12 
DO 15 K=l,100 
STAR=1 

RO ( 3 ) = (R0( 1 )--DVOL( 2 ) ) / ( DVOL ( 2 ) ~DVOL (1 ) ) + ( RO ( 2 ) *D VOL ( 1 ) ) / (DVOL( 1 ) 
“DVOL( 2 ) ) 



GO TO 01 

I F ( ABS( DVOL { 3 ) )“1» F“04 ) 17*17» 1 
I F (DVOL ( 3 ) *OVOL (2) )14»17»15 
RO( 1 )=RO(3 i 
DVOL( l)=DVOL( 3) 
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appei^:dix a 



Input Data 

In orgenizin" the machine program that was ultimately 
used, it x-ias necessary to anticipate what information v;as 
likely to be available or fixed for the machine design. It 
was anticipated that certainly the desired power output 
would be iQiovm. The angular velocity of the rotor and 
povrer factor of the load may also be determined from 
particulars of the intended machine application. Further, 
it was expected that materials selection would precede this 
portion of a design, therefore various material character- 
istics x-xill be available. Also expected is some knovrledgo 
of construction techniques to be employed, e.g. the type of 
armature cooling used and characteristics of the rotor 
design. Finally, some information concernixig the dewar is 
necessary. Specifically the anticipated thickness of the 
deT-:ar vrall must bo provided. Thus, input data will include: 

a. pov;er output, P^ 

b. mechanical frequency, Wm 

c. maximum armature current density, 

d 

d. dovrar wall thickness, A 

e. maximum tip spood of dewar, 

(fixed by strength considerations in rotor design or 
limitations imposed by dewar technology) 

f. po\i7er factvor, pf 



II 



t 
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g. flti.x dens5.ty to cause shield saturation, 

h. average density of rotor, armature and shield, p^, p^ 

As mentioned in the Introduction, the sujpercondvicting 
material used in the field windings will have characteristic 
values of critical current and critical fliox density. There- 
fore, additional inuut data will be a value of ciu?rent density 
in the field windings, J^, and the associated value of critical 
flux density, B 

max 

For any field winding design, the geometry of the magnetic 
fields can be found. If is knovm, the maximum flux density 
can bo calculated. This maximmi flux density must be found 
to cxisurc that it does not exceed the critical flux density 

associated v:ith (B ). 

.(K a X 

Convci'sely, if and B jti ^ are given, certain 
limitations on field vjinding design are imposed. 

It v;as found expedient tc provide a body of data irhich 
could be vised in the computer program for finding the best 
field winding design fox* some given and B The 

method for obta.ining and utilizing this data 
in Appendix B. 



is described 
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APPEIIDIX B 



Input Data for Field V/indinp; Design 

In the design of the field v;indings, three assumptions 
are made. First it is assumod that the fundamental component 
of the magnetic fields is the significant part. This assump- 
tion is fairly standard and has been shovm to be acceptable . ^ 
Second, it is assumed that the fields generated by the arma- 
ture may be neglected. Since the armature fields oppose the 
fields generated by the field windings, this approximation 
is conservative. Third, it is assximed that the fields 
reflected by the shield may be neglected. Althougli these 
fields add to those from the field vrinding, they are small, 
becoming insignificant as p increases. These approximations 
are made in order to deer ease the nimiber of variables in the 
calculations . 

From Table I, it is seen that in the region Z p Z 
the maximum field intensity will be in the p -direction 
for 0=0. In accordance v;ith the above approximations, 
for p r 2, 
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e " 






“2p {2 • 


• P) 


^2 /RJ 

'■“I - 2-pj 


p+2 






lei i e 1 





(B-1) 



or 



B 

const . 



Bptr 







+(2-p) 




p+2 

-(2-p)y X 




p-i 1 ■' 



0 / 



(B--2) 
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I, 
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I 



\ 



If, for a given value of y. 



-k2- 



d(B p /const . ) 

dTpTRp 



= 0 



(B-3) 



is solved for ( p /R^) , and the value obtained used to 
evaluate Eq. (B-2), the iriaximum value of Bp /const, is 
obtained for that y. 

For the case p = 2 
2J 

In • 

9 



H = 



TT 



1 '*1 

in 



(B-l,.) 



and 



B p rr 



const . 



-/*o'^f^2 



1 £ V 

IT y 



3\ 



P 1 

iij, e' 



Thus, if is known, it is possible to obtain 



(B-5) 



B 



P max , 



f(p,y) 



(B-6) 



In the follov;ing computer program, for the IBil 1130, 

B /const, is computed for values of o = 2 bo 6 and 

^ ruQ--^ 

y = 0 . 05 ^ to 0.95* I'tio results arc displayed graphically 
in Figoiro 10. Using these results, it is possible to find 
the value of y corresponding to a given B . /const., 
thvis fixing the radial dimensions of the field vrindings. 
The case p - 1 is not computed since for this value of p 



B 



mac 



/cor:, t, - 1 - 



(B-7) 



In the program, the follovring notation is used: 



p = number of pole pairs 



y = 

R = 9 /^ 2 ^ 
of p /R 



is the inside field radius, etc. 

i.e. R1 in the program is the first value 
2 and is not the inside field radius. 



B = Bp /const. 



B MAX = B ^ / const, 

p max 



The program procedure is simply a solution of Eq. (B~3) 

and evaluation of Eq, (B-2) or Eq. (B-5) fo^ the ma,ximum 

( 7 ) 

value. For the case p « 2 a false position method ' is 

used to solve Eq. (3-3) > whereas the Kevrton-Raphson 
( R ) 

techniqiie' ' is employed for all cases p 2. 

The printed results of this program include values of 
p, y and R p v'hich aro used in the design program. 

The data reprosonus the curves in Figure 10, Within the 
design program, P3 ^ /const, is computed, vrith vzhich it is 

P TflrtJC 



then possible to obtain from the curves a value of y ’A'hich 
will satisfy the constraint imposed by the cairrent and 
critical flux density. 



The program -used to solve for the data in the array 
R(p>y)> subrou.tine used to interpret B(p,y) to 

find a val’ae for y, are incli.uded in the foil oaring pages. 
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BMAX = , 25*RNEW - « 2 5^ ( Y^'^4 ) / ( RNEW**3 ) ~ RNE W^ALOG { RNEW ) 

GO TO 90 

USE NEWTON-RAPHSON METHOD TO FIND RADIUS AT MAX RHO DIRECTED B FOR 
CASE P NOT EQUAL 2. 



A= ( P->:-*2 + P-2 )* (R**( 2*P) ) ( R** ( P + 2 > )- < P**2-P-2 ) * ( Y** { P+2 ) ) 
B^2-:! P-“‘- ( P*->2+P“2 ) ( 2 '*-P-l ) ) -2*P--- ( P + 2 ) ( R’f-* ( P + 1 ) ) 

RNEW=R - A/B 
CHECK CONVERGENCE 

IF(ABS(R - RNEW) - l.E-04) 23*23, 21 
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DATA READ TO THE MACHINE AS ARRAY 
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APPENDIX C 



Pesi|5;n Pro;!;raia 

The design computer program, i^rritton for the IBM 1130, 
is an iterative solution of the basic design equations to 
find that machine design which requires the least volume. 

The program uses the normalized variables x and y which 
represent the thickness of the armature and field vrindings. 
The variables and R serve to fix the radial location 



of the windings. The machine length is found from consider- 
ation of the working volume required to obtain the necessary 
internally generated pov;er. The reductd.on from internal 
power to rated povjor is obtained from Eq. (11-15). The 
iterative procedure followed is based on varying values 
of R , which vras fovmd to be the most convenient parameter 



2 



o 



O' 

fo'r this purpose 



Since 



Vol. = TT R^^ 1. 

OS t 



(C-1) 



the variation of R 



w'ith R^ was first examined. By 



os 



evaluating Eq 



obtains 




. (xI-21) and substituting in Eq 

/w \ P'2 



(11-22) one 



R 



os 



o 



(C-2) 




R 



o 



(C-3) 




! 
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Flgiipe 11 Vai'iation of vzith 

For small values of R , the shield is in the vicinity 

o ^ 

of intense magnotic fields and R must be largo. As R 
increases, the thickness of the shield decreases but R 

os 

is still larger than R^. The value of R^ v;hich requires 

the smallest R is found to be 
os 



R 



o opt . 



(pc^) 



1/p+l 



(c-4) 



Though this value of R^ does not necessarily yield a 
minimum volume machine (since no variation of 1. is 
considered), it is generally near the rnininum volume value 
because of the significance of R ^ in Eq. (G-l). 

An obtained from Eq. (C-if.) vrill satisfy the constraint 
of Eq. (11-22), However, it must also satisfy Eq. (11-20). 
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Since 



= iTT??) 



u.^«> (M 



p-^2 



R ^ B . , 

o — rated 



(c-5) 



it is required that 



Ro ^ (C^/p) 



1/p-l 



(C-6) 



Since 



opt. “ “ p * 



> .-1 



(c-7) 



Eq. (J.I-20) is also satisfied. 

Turning now to the design program, the follovring 
procedure is carried out for successive values of p: 
From the input da.ta., one obtains 

R. = V. . /o) , 

1 txp ^m’ 



(C-8) 



Rp = - A . 



(c-9) 



Evaluating 
B 



B , , rr 

/const, r 

ma.>.' ‘ ‘ " 2 u "’J/.H'T 

^ o X 2 



(C-10) 



and entering the curves of B /const, vs. y, described 

in Appendix B, a value of y is obtained, -vith this 

information it is D0ss5.ble to evaluate C., and find R ^ 

^ 1 o op o . 

vihich is used to ca.lculate an initial machine design. 

Having selected a value for R one obtains 
^ o 



X 



R. 

X 

j-f" 



(C-11) 










i 




At this point in the program, the differonco bctvreen tho 
internally generated pov;er and the desired pov;er output is 
unlmovm. Assuming that this difference is largely due to 
load power factor, it is approximated that 



= P^/pf 
g o 



(C-12) 



Since 



’g =iVa 



(c-13) 

(c- 14 ) 



Inhere 



CO„ = pCO 



m 



Ti (l-y^)Pv| 



(C-l5) 

(C-16) 



a 



6N 



(C-l?) 



and the appropriate expression for M is selected from 
Table II, one obtains 



S 



o'-'^m-af (l-y^*^) 



TT 






p-^2 






l~x“P -2 






k~r/ 



(2'P)' 



R 



o 



for p ^ 2, or 



(0-18) 
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P =4 

g 2 



(1-y^) r| “ 






1-x' 

T-” 



4 



(C-19) 



for p = 2. P^’ora Eq. 

solve for 1 . Then, 
m 



Eq. (11-19), 



1 



t 




= 1 



-f- 

rti 



tr 

3p 



(C-l8) or Eq. 
in accordance 



(l+x)R^. 



(0-19), it is possible to 
with Eq, (II-18) and 



(C-20) 

(C-21) 



It is noxT possible to find the reduction in power fr-oro 

P to P . Applying the law of cosines to tho diagram in 
6 ^ 

Figiu’e 5 results in 






a 






Normalizing, one obtains 



(V /V -5 
t a a 



2 V, /V X coJ 
t a a 



(| 6) 



-1 = 0 



(C-23) 



v;here 



x^ = x^ I /V (C-24) 

a s a a % -r/ 

is the per unit araiature reactance loss. If Eq. (C-23) 
is solved for V^/y^, then Eq. (II-15) can be evaluated 



for P . 
o 

Tlio value of P thus obtained 
o 

value specified in tho input data, 
grvsatsr than a given tolerance (1 ®j 
for this urogram), the va.luo of F 



is compared with the 
If the difference is 
of tho desired value, 
assvned in Eq, (C-12) 
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is increased (oi* decreased) by the percentage difference 
betvreen actual and desired oiitput power (see statement 30 
in the design program). The program then returns to Sq. 
(C-l8) or Sq. (C-19)> and calculation continues in this 
fashion \mtil an acceptable difference is reached. 

VJhen a satisfactory is obtained the program contimies 
with calculation of 




o 



(C-3) 



Vol. = rr 1. 

os t 



(C-1) 



VJeight = TTl^( (l~y^)R 2 p 



f 





(G-25) 





Tr^(l-x^)Rp 

o 





for p 5 ^ 2, or 




X 



(C-^27) 




i 
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for p = 2, and 



N I = 
a a 



Tr(l-x2)R2 



(C-28) 



Tho calculations described above are carried out for the 
initial value of R^. The program continues by incrementing 
Rq (by 0,1 in this program), returning to Eq, (C-11) and 
repeating all the steps through Eq. (0-28). R^ continues 
to be incremented until it is found that the volume for 
successive designs is increasing, which indicates tha.t the 
minimvim voluitie design is within the range of those designs 
that have been completed. The program then calls a svi.b- 
routine, described in Appendix D, Tdiich evaluates tho 
data for the designs that have beon completed, and from 
thifi data determines the mini mum volu^aO des.Lgn. 

V'ithin the design program, the following notation is used 

PV.’R = desired P 

o 

POUT = actual P 

o 

B(p»y) - array of data representing the graphs 
in Figure 10 

OUTPT K array in vjhich parameters of each design 
are stored 



BC = B 



p max/const , 



XA = un-normalised amat’are reactance, v;hich 
becomes normai.iaed when stored in OUTPT. 



VAOII = V /N 
a a 



r / >T 

V a 



VTCH 
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PG = Pg 
ROS = R 

os 

VOL “ Volume 

V/GT = \foight 

HIA = N I 
a a 

A print-out of the souj?co deck for the design program 
is included in the follov.'ing pages. 
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APPEKDIX D 



Subroutine to Find the Mininnim Voliitne Design 

After the design program computes throe designs of 
successively increasing volume, it calls the subroutine 
MIN, The data stored in OUTPT at this time will include 
at least tv;o designs beyond the minimum vo3ume design. 



Vol. 



^opt ^01 ^02. ^03 ^04 

Figure 1 ,2 a. Data Stored in OUTPT 
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The program computes dVol/dE^ for the three values of 
that bracket the mini.mum volume design. 

Since 



Vol = tf RX^ 1. 

os t 



(K-1) 



then 



dVol 

Hr 



= rr R 



2 

os dR 



+ tT"! — — 

^•^t dR 



d(R^ ) 
os 



(D-2) 



Considering Equations (C-l8), (C-19)> and (C-21), it is 
possible to express 



1. = ^ rP"^ + (R R. ) 
t o ,3 p o 1 



(D-3) 



\jhere 



f(x). 



(D-4) 



Therefore 



dl. P o P dC ^ 

= (P“2) rP"^ § rP"^ 

dj^ M' / Q ^ „2 dR o 

O O 0 o 

c 



-f 



T( 



3p 



(D-5) 



and, from Eq. (C-3), 



os • 

dR ““ 



-2pcfR’^P"^ » 2C- (i-t))rrP 

^ 1 o 1*0 



2R. 



o 



(D-6) 



When values of DVOL for R^^ , have been 

computed, linear interpolation is iised to find the value 



of R for xdiich DVOL (R ) 0. Control then returns to 

o o 

the main design progr.fim xdi.ero the design characteristics 
for this v^;.lue of H 



arc found . 
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A print-out of tho soui-ce dock for this sub-routine is 



included in the 



f outgoing pages. 
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APPEimiX E 



Effect of the Shield on I-Iaximxjjn Flux Density 

As described in Appendix B, the assumption was madCj 
in determining characteristics of the field vrindings, that 
the effect of the shielding could be neglected. The 
significance of this assumption is examined as f ollovxs : 

If B p /const, is computed, as in A.ppendix B, vrithout 
neglecting shield effect, the result is 



B, 



Const . 






-2p 



R. 



(2-rp) 



p\P-l 

R 



-(2-p)yP"2^ 



R. 

[T 



p^i 



(2-p)(l-y^‘^'^) 



for p / 2, en.d 



B 



/ 



Const . 



1 

If 



P 4 

it: “ ^ 






M*'' p 



Recognizing that 

R. 
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«i - A 



R 
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vjhero 
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r“ 
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p-1- 







IP 
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- 8 



(S-l) 



R 



In 



2 ... 






(S-2) 



(S-3) 



(s-4) 



Eq, (E-1) and Eq. (E-2) can be vrritben in tox.-s of th: 
nor.mali zed nara.’ ioters x, *r and cS . 
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Assxuning that S is constant over a range of machine 

sizes, it is possible to compute the actual B ^ /const. 

> ^ ^ p max 

in terms of x and y. The solution technique is the same 

as that described in Appendix B. One can then compare 

the value of B /const, obtained from the simplified 

p max 

equations with that value obtained from the method 
described above. 

The results of such a comparison are displayed in the 
following graphs. The data is for the cases of p = 3 to 6, 
which are of most interest. 
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Figure 13 
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Figuro 16 
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